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Abstract 

Measurements  of  the  wave  speed  c  in  Ml  imply  a  fast  lime  A  =  it-/ pc2  of 
relaxation.  This  and  the  delayed  die-swell  measurements  suggest  that  Ml  is 
not  very  elastic.  Extensive  and  very  reliable  values  of  the  climbing  constants 
show  that  Ml  has  weak  normal  stresses  at  the  level  of  STP.  The  climbing 
constants  plus  back-extrapolated  values  of  jV,/y2  taken  from  measurements 
at  two  temperatures  by  Binding,  and  Walters  and  Prud’homme  lead  to 
values  of  the  two  coefficients  in  the  second-order  regime  of  slow  flow  of  Ml. 
This  gives  us  all  the  constants  in  the  Roscoe  formula  for  the  quadratic 
correction  of  Trouton’s  viscosit}'  and  also  allows  us  to  compute  the  normal 
stress  ratio  —  A^/JV,  =  0.11,  independent  of  temperature. 

Keywords:  climbing  constant;  delayed  die  swell;  cxlcnsional  viscosity;  wave  speed;  test  fluid  Ml 


1.  Introduction 

The  rheometrical  data  presented  in  this  report  are  unique  in  that  it  is  the 
only  report  on  the  properties  of  the  Ml  fluid  to  present  measured  values  of 
the  climbing  constant,  of  the  two  constants  of  the  second-order  approxima¬ 
tion  of  a  viscoelastic  fluid,  of  the  intensity  factor  for  extensional  viscosity, 
and  secondary  motions  at  second  order.  We  computed  N2/Nx  =  —0.11  at 
low  shears  from  our  measured  climbing  constant  at  20  °C  and  back  extrapo¬ 
lation  of  iV,/2y2  given  by  K.  Walters.  We  computed  N2/Nx  =  —0.12  at  low 
shears  from  our  measured  climbing  constant  at  27.2°  C  and  back  extrapola¬ 
tion  of  Nx/ 2y2  as  given  by  R.K.  Prud’homme.  We  are  also  the  only 
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laboratory  to  give  values  of  the  wave  speed,  the  rigidity  and  the  effective 
time  of  relaxation.  Our  time  of  relaxation  is  shorter  than  that  computed  in 
the  traditional  way  using  standard  rheometers.  The  reason  for  this  is  that  the 
response  time  of  the  wave  speed  meter  is  much  shorter  (say  10~4-10-3  s) 
than  the  response  time  of  cone-and-plate  rheometers  (10-1  s).  The  speeds 
which  we  measure  and  the  associated  relaxation  times  seems  to  correlate 
with  various  critical  phenomena  which  occur  in  the  flow  of  polymeric 
liquids,  the  most  notable  being  delayed  die  swell.  We  present  some  data  for 
delayed  die  swell  in  the  Ml  solution.  The  Ml  solution  has  a  fast  time  of 
relaxation  as  is  characteristic  of  fluids  with  a  rather  more  Newtonian  than 
elastic  response. 

2.  Theory 


The  operational  theory  for  the  measurement  of  the  climbing  constant  is 
given  in  the  1984  paper  of  Joseph  et  al.  [1].  The  height  function  for  the  climb 
is  given  by 

H(r.  Q)  =  h0(r)  +  h2{r)&,  (2.1) 

where  £2  is  the  angular  frequency  and  terms  of  0(&4)  have  been  neglected. 
Here  h0(r)  is  the  static  climb,  independent  of  8,  and  h2{r)  is  to  be 
determined  by  the  following  problem: 


4  4 

a  n  .  P° 


-(rhi)  -psK=  -2-^-/?  +  -^-, 

h2(a )  =0,  h2(r )  -»0  as  r-*  oo 


where  T  is  the  surface  tension, 

fi  =  3a,  +  2  a2  (2.3) 

is  the  climbing  constant,  a,  and  a2  are  constants  of  the  quadratic  approxi¬ 
mation  to  S  in  the  stress  T=p\  +  S,  S  =  5,  +  S2 . . . , 

Sx=  \lA-\U\  is  the  linear  or  Newtonian  approximation, 

Si=<X\A2[U]  +<x2A][U).  (2.4) 

The  expression  S(2)  ~Sl  +S2  is  called  the  extra  stress  of  a  fluid  of  grade 
two.  This  is  just  a  way  of  speaking;  there  is  no  fluid  of  grade  two,  but  such 
an  expression  arises  on  every  simple  fluid  when  the  flow  is  slow  and  slowly 
varying;  for  a  recent  full  discussion  see  Joseph  [2].  The  constants  a,  and  a2 
are  related  to  the  zero  shear  values  of  the  first  and  second  normal  stress 
differences  and  N2  by  the  following: 


r„-7’22  =  Wl  =  -2a,y2, 

T22  —  r33  =  N2  =  (2<X|  +  a2)  y  , 


(2.5) 


where  y  is  the  shear  rate  and 

7l2  =  f/t- 
If  we  define 

«,=  limiVj/y2, 

f-*0 

«,  =  IimjV-j/y2. 

f-0 

then 

«i  = 

a2  =  it1  +  n2. 

/?  =  3a,  t  2a,  =  -j-  +  2n,. 


(2.6) 

(2.7) 


(2.8) 


The  fluid  will  not  climb  a  rod  if  n2/nx  <  —1/4.  The  theory  of  rod 
climbing  with  surface  tension  neglected  was  done  by  Kaye  [3],  Joseph  and 
Fosdick  [4],  Hoffman  and  Gottenberg  [5],  and  Bohme  [6].  Joseph  et  a!.  [7] 
were  the  first  to  use  theory  in  experiments  designed  to  measure  fi.  They 
showed  that,  to  get  the  predicted  shapes  of  the  free  surface  to  agree  with 
measured  ones,  it  is  necessary  to  retain  the  effects  of  surface  tension,  as  in 
(2.2).  A  very  accurate  approximate  solution  of  (2.2)  was  derived  by  Joseph  et 
al.  [7j.  When  evaluated  at  the  rod  this  solution  gives 


/  (  \  a  (  W 

hi0)  2TyfS  (  4  + A 


pa 2  \ 

2+\y 


(2.9) 


where  A2  =  a2S  and  S  =  pg/T.  Joseph  et  al.  [7]  and  Beavers  and  Joseph  [8] 
proposed  (2.2)  and  (2.9)  as  a  basis  for  a  rotating  rod  viscometer.  A  possible 
defect  of  this  kind  of  rheometer  is  that  it  requires  values  for  the  surface 
tension.  Fortunately,  the  values  of  T  which  have  been  measured  in  the 
different  liquids  used  in  experiments  and  in  Ml  are  nearly  the  same,  about 
30  X  10-5  N  cm-1.  An  error  analysis  given  by  Joseph  et  al.  [1]  shows  that 
the  value  computed  for  fi  from  the  graph  of  the  measured  height  of  climb  at 
the  rod  is  not  strongly  affected  by  small  changes  in  surface  tension. 

Roscoe  [9]  presented  a  formula  for  the  extensional  viscosity  at  low  rates 
of  shearing 

Tn  -  T22  =  3[sfi+  (a,  +  oc2)s2]  +  0(\s\3),  (2.10) 

where  s  is  the  rate  of  stretching.  This  is  probably  the  only  model-indepen¬ 
dent  formula  for  extensional  viscosity  which  can  be  evaluated  numerically 
by  relatively  simple  rheological  measurements.  Equation  (2.10)  was  given,  in 
a  slightly  different  form,  by  Zahorski  [10]. 
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Using  (2.8).  \vc  find  that 

a\  "**  az  ~  ~T  +  n2-  (2.11) 

Fluids  with  n2/nl  <  — 1/2  stretch  weaken.  Equation  (2.8)  3  show's  that  such 
a  fluid  will  not  climb  a  rod.  We  also  note  that 


a.  =  ln2  -  /?. 

<r,  =  2/S-3»,,  (2.12) 

al  ***  a2  =  P  ~  11 2- 


Equation  (2.12)3  shows  that  the  intensity  factor  for  extensional  flow',  the 
second-order  correction  of  the  Trouton  viscosity,  is  equal  to  the  climbing 
constant  when  n2  =  0.  It  is  easily  verified  that  a,  +  a,  —  jfi  when  n-Jnx  = 
-1/10. 

Joseph  [2]  has  shown  that  the  intensity  factor  a,  +  a,  for  the  extensional 
viscosity  at  second  order  also  controls  the  creation  of  vortices  at  second 
order  in  the  following  sense:  suppose  that  the  vorticily  £  —  +  e2£2  +  0(e3), 

where  €  is  a  parameter  perturbing  rest.  Then  we  find  that  £,  and  £2  satisfy 

V2£j  =  0  (2.13) 

and 


3$, 

p~d 7 


+  p  curl(£,  A  ux)  —  p  v2£2  +  +  tf2)cur^  div/12. 


(2.14) 


Equations  (2.13)  and  (2.14)  show  that  the  effects  of  viscoelasticity  on  the 
vorticity  appear  first  at  second  order  with  an  intensity  factor  a,  +  a,. 

How  should  the  constants  a,  and  a2  be  measured?  The  climbing  constant 
(2.3)  gives  one  combination  of  cq  and  a2  which  can  be  measured  with  very 
good  precision.  The  hard-to-measure  second  normal  stress  is  related  to  cq 
and  a2  by 

n2  —  2a, +  a2.  (2.15) 

It  is  claimed  that  n2  can  be  measured  by  measuring  the  surface  deflection  in 
a  tilted  trough.  A  full  discussion  of  this  method  is  given  by  Tanner  [11],  The 
shear  rates  which  can  be  achieved  in  a  trough  tilted  through  a  small  angle 
are  rather  small  so  that  the  deflections  which  develop  on  the  free  surface 
would  not  be  easy  to  measure  in  many  fluids.  The  comparison  of  Table  2 
and  Table  3  giving  cq  +  a2  for  the  fluids  D1  and  Ml  suggests  that  the 
surface  deflections  in  Ml  would  be  more  than  50  times  smaller  than  in  Dl. 
In  fact,  rheological  data  using  this  method  are  rather  sparse. 

The  second  normal  stress  can  be  obtained  in  measurements  on  ordinary 
cone-and-plate  rheometers  by  a  well-known  method.  The  problem  with  this 
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method  is  that  the  values  of  the  first  and  second  normal  stresses  at  low 
shearing  cannot  usually  be  obtained  because  the  transducers  will  not  work 
unless  the  shear  is  above  a  certain  threshold.  However,  it  is  very  easy  to 
show  [2]  that,  at  low  rates  of  shearing,  the  ratio  of  the  normal  stress 


differences  has  a  limiting  value 


A',  1  ’ 


a2 

2<r, 


+  0(f), 


(2.16) 


with 


1  + 


J-a, 


(2.17) 


Hence,  we  may  hope  that  (2.15)  may  be  measured  by  backward  extrapola¬ 
tion  of  the  values  of  —  JV-j/jV,  measured  on  cone-and-plate  rheometers.  The 
values  of  and  a2  may  now  be  determined: 


1 

1  4-  4/2  2/>;, 


P, 


(2.18) 


„  1  +  n2/nx  - 
“2  2l  +  4  n2/nl^ 

1  +  2 h2//2,  . 
“'  +  a^i  +  4 


(2.19) 

(2.20) 


Evidently  N2/2y2  is  too  small  to  measure  in  Ml  so  that  back  extrapolation 
of  N2/Nx  without  considerable  scatter  cannot  be  carried  out.  However,  it 
seems  that  the  graphs  of  Nx/2y 2  shown  in  Figs.  8  and  9  can  be  extrapolated 
to  give  —a,.  Then  a,  and  ft  determine  a2  and  N2/Nx  uniquely. 


3.  Climbing  constants 


The  climbing  constants  (Fig.  1)  were  measured  using  a  rotating  rod 
viscometer.  We  have  figures  like  those  of  Figs.  2-5  for  every  one  of  the  38 


Fig.  1.  Summary  of  the  dependence  of  climbing  constant  4  on  temperature  T  (6.2-41 .5  °C). 
The  equation  of  the  line  is  =  39.7  expf-O.lSSr)  g  cm-1.  At  T=  20 °C,  /?  =1.68  g  cm-1. 


0 


0.5  0.6 


0.1  0.2  0.3  0.4 

o2  (rev/s)2 

Fig.  2.  Height  rise  h  (cm)  at  the  rod  vs.  the  angular  velocity  squared  or  ((rev  s-1)2).  o 
(angular  velocity  increases).  h=  —0.011  +  0.817&r;  □  (angular  velocity  decreases).  /z=1.8x 
10-3  +0.813«2.  The  average  climbing  constant  is  ft  =  23.2  gem-1  at  T  —  6.2° C. 


Fig.  3.  Height  rise  h  (cm)  at  the  rod  vs.  the  angular  velocity  squared  or  ((rev  s-1)2).  o 
(angular  velocity  increases),  It  =  6.65  X 10“ 3  +0.0478 to2;  □  (angular  velocity  decreases),  h  — 
6.99x10~3+0.0479w2.  The  average  climbing  constant  is  p  =  1.43  g  cm-'  at  T =19.45°C. 


Fig.  4.  Height  rise  h  (cm)  at  the  rod  vs.  the  angular  velocity  squared  w2  ((rev  s-1)2).  o 
(angular  velocity  increases),  h  —  9.00X10“ 3  +  0.003  35 w2;  □  (angular  velocity  decreases), 
h  =  7.70X10-3  +  0.00338w2.  The  average  climbing  constant  is  /?  =  0.168  g  cm'1  at  T= 
35.55°  C. 
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a 2  (rev/s) 2 


Fig.  5.  Height  rise  h  (cm)  at  the  rod  vs.  the  angular  velocity  squared  ur  ((rev  s-1)2).  o 
(angular  velocity  increases),  for  the  first  four  points  /»  =  1.19  Xl0-3  + 7.90 Xl0_4w2;  □ 
(angular  velocity  decreases),  for  the  first  four  points  /;  =  0.0133 +6.53x1 0-' ’w2.  At  this 
temperature,  the  fluid  sinks  for  higher  values  of  angular  velocity.  The  average  climbing 
constant  is  /?  =  0.093  g  cm-1  at  T  =  41.5°C. 

entries  in  Table  1.  In  Table  1  we  also  list  the  test  date  as  a  check  to  see 
whether  the  Ml  is  time  dependent  or  not.  As  indicated  in  the  table  Ml  is 
not  time  dependent,  at  least  for  a  period  of  time. 

Surface  tension  of  the  fluid  T  =  29.9  X  10“5  N  cm-1,  density  of  the  fluid 
p  =  0.895  g  cm-3  and  the  radius  of  the  steel  rod  a  =  0.474  cm.  T  and  p  are 
almost  constant  for  the  temperature  range  of  our  measurements.  The  surface 
tension  was  measured  with  a  ring  tensiometer. 

The  diameter  of  the  fluid  container,  which  is  made  of  copper,  is  10  cm. 
The  angular  speed  of  the  rod  can  be  measured  to  within  0.002  rev  s-1.  The 
measurements  of  the  height  of  climb  are  repeatable  to  within  0.002  cm.  The 
rod  and  fluid  are  enclosed  in  a  Plexiglas  chamber  which  is  connected  to  a 
temperature  controller.  The  temperature  is  measured  at  the  surface  of  the 
fluid  and  the  probe  is  kept  at  the  same  position  for  all  the  measurements. 
There  are  temperature  differences  between  the  fluid  and  air  in  the  chamber; 
for  more  accurate  measurement  we  need  a  better  way  to  control  the 
temperature. 

At  larger  rotational  speeds  the  steady  climbing  bubble  loses  stability  to  a 
time-dependent  motion  which  has  been  called  the  “breathing  instability”. 
This  motion  first  appears  at  a  critical  w  =  coc,  then  grows  quickly  to 
time-periodic  motion  with  a  long  period  in  which  a  band  on  fluid  rises 
slowly  to  a  certain  height,  then  collapses  down  into  the  body  of  fluid.  The 
time-periodic  breathing  motion  in  Ml  is  shown  in  Figs.  6(c)-6(f)  and  in 
photographs  displayed  in  Figs.  95.2  and  95.3  by  Joseph  [12].  A  general 
description  of  the  breathing  instability  is  in  the  reference  just  mentioned,  in 
Beavers  and  Joseph  [8],  and  in  Joseph  and  Beavers  [13]. 

The  first  quantitative  data  for  the  breathing  instability  are  presented  in 
Fig.  7  where  the  critical  speed  is  plotted  against  the  temperature  and  fit  to 
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TABLE  1 


Measured  values  of  climbing  constants.  AH  data  was  taken  in  February  1989.  Ml  is  not  time 
dependent  at  least  for  a  period  of  time 


Temperature 

(°C) 

Climbing  constant 

P  (gcm_!) 

Date 

6.2 

23.2 

14 

6.3 

21.8 

14 

9.15 

11.3 

15 

9.15 

11.5 

15 

12.25 

5.80 

15 

12.3 

5.76 

15 

16.3 

2.80 

16 

16.3 

2.87 

16 

18.9 

1.67 

16 

18.95 

1.59 

16 

22.1 

0.931 

16 

22.2 

0.934 

16 

25.4 

0.568 

17 

25.3 

0.585 

17 

29.05 

0.366 

18 

28.95 

0.379 

18 

32.9 

0.219 

18 

32.95 

0.217 

18 

37.25 

0.130 

19 

37.35 

0.128 

19 

41.5 

0.083 

20 

41.5 

0.093 

20 

35.5 

0.173 

20 

35.55 

0.168 

20 

31.15 

0.304 

21 

31.15 

0.317 

21 

26.9 

0.576 

21 

26.95 

0.564 

21 

23.15 

0.876 

21 

23.1 

0.873 

21 

19.45 

1.43 

22 

19.45 

1.43 

22 

18 

1.97 

22 

17.95 

1.88 

22 

13.9 

4.61 

23 

13.8 

4.53 

23 

11.45 

8.49 

23 

11.35 

8.75 

23 

the  curve  wc  =  exp(0.074!T).  Recalling  now  that  /?  =  39.7  exp(  —  0.1587)  we 

I  A 

may  conclude  that  the  critical  speed  wc  is  proportional  to  l/\//?,  approxi¬ 
mately. 
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4.  Limiting  values  of  the  normal  stresses 

Thrust  measurements  on  the  plate  of  a  cone-and-plate  rheometer  give  the 
normal  stress  difference  TV,( y)  as  a  function  of  the  shc«;  rate  y.  The  same 
thrust  measurements  in  a  parallel  plate  rheometer  give  values  of  N2  -  TV, 
where  TV2(y)  is  the  second  normal  stress  difference.  Neither  measurement 
can  be  easily  carried  to  low  rates  of  shear  because  of  the  insensitivity  of  the 


Fig.  6(a-f)-  Photographs  of  rod  climbing  in  Ml  at  21.5  °C.  (a)  At  rest,  there  is  a  small  static 
dim'  due  to  capillarity;  (b)  stable  steady  climb  co  =  3  rev  s_1;  (c)  and  (d),  slightly  post 
critical  time-periodic  motion  at  <o  =  4.94  rev  s_1;  (c)  and  (0,  large  amplitude  time-periodic 
motion  at  <o  =  6.55  rev  s_1. 


transducer.  From  these  two  measurements,  we  can  compute  N2/N1  as  a 
function  of  y  and  try  to  extrapolate  it  to  zero  shear.  This  first  procedure 
works  well  for  the  test  fluid  D1  used  in  the  second  normal  stress  difference 
project  [14].  From  Fig.  5  of  Ref.  14  we  easily  estimate  by  extrapolation  to 
zero  shear 
N 

0.07  <  -  -^  <  0.10  at  25  °  C.  (4.1) 

The  climbing  constant  for  D1  at  25  °  C  was  measured  by  us  as 
0—  62.2  g  cm-1. 


(4.2) 
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In  Table  2  we  have  given  the  values  of  the  quadratic  constants  cq  and  a2 
computed  from  (4.1)  and  (4.2)  using  (2.8)  and  (2.17)  and  the  values  of  the 
limiting  normal  stress  /q  and  n2  using  cq  and  a2  in  (2.5)  and  (2.7). 

The  normal  stresses  in  Ml  are  more  than  50  times  smaller  than  in  Dl. 
However,  it  appears  to  be  possible  to  back  extrapolate  Nx/ 2y2  to  its 
limiting  value  —  cq.  This,  together  with  the  climbing  constant  =  3a,  +  2a2, 
gives  a,  and  a2  separately.  Figure  8  gives  measured  values  of  Nx/2y2  for 
Ml  at  20  °C  which  were  measured  in  Professor  Walters’  laboratory  on  a 
Weissenberg  rheogoniometer.  Figure  9  gives  measured  values  of  Nx/2y2  for 
Ml  at  21.2°  C  which  were  measured  on  a  System  IV  Rheometrics  rheometer 
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Tea-re? 

Fig.  7.  Dependence  of  the  critical  speed  of  instability  uc  on  temperature  (°C).  The  equation 
of  the  line  is  ttc  =  exp{0.0747')  rev  s  *. 


TABLE  2 


Material  parameters  for  D1  at  25° C  using  (4.1)  and  (4J2) 


—  =  0.10 

—  »,/«,  —  0.07 

/?(g  cm-1) 

62 2. 

611 

a,  (gem-1) 

-103 

— S6 

a,  (gem-1) 

1S7 

161 

a,  +  a. 

84 

75 

//,  (gem-1) 

206 

172 

//,  (gem-1) 

-17 

-11 

with  a  sensitive  transducer  in  the  laboratory  of  R.K.  Prud’liomme.  From  the 
extrapolated  values  a,  =  —3.0  g  cm-1  at  20°C  and  a,  =  —1.0  g  cm-1  at 
27.2 ° C  and  the  values  0  =  1.68  g  cm-1  at  20 ° C  and  0  =  0.54  g  cm “ 1  taken 
from  Fig.  1,  we  compute  the  entries  in  Table  3. 


Fig.  8.  Measured  values  of  Ny/ 2y2  for  Ml  at  20  °C.  The  backward  extrapolation  gives  that 
N^/ly2-*  -  «,  is  about  0.3  Pa  s2  or  3  g  cm-1.  (After  Binding  ct  al.  [15].) 
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Fig.  9.  Measured  values  of  ;V,/2y:  for  Ml  at  27.2°  C.  The  measured  values  of  the  stress  for  O 
are  below  ihe  allowed  transducer  limits  and  are  discarded.  We  took  the  limiting  value  to  be 
—  or,  =  0.1  Pa  s2  or  1  g  cm  '  *.  (After  Prud*homme  (16{.) 


Table  3  shows  that  As/Af,  is  very  insensitive  to  changes  in  the  tempera¬ 
ture  even  over  a  range  of  temperature  in  which  the  values  .V,  and  jV,  depend 
strongly  on  the  temperature.  The  value  of  As/iV,  is  also  not  very  sensitive  to 
errors  in  the  measurement  of  «,  =  —2a,.  At  the  Combloux  meeting  mea¬ 
surements  at  20  °C  were  reported  in  the  range 

2.77  g  cm ' 1  <  —a,  <  3.50  g  cm'1. 

Using  P  =  1.68  we  compute 

0.097  <  -NJN,  <  0.129. 

This  range  of  values  is  consistent  with  rV,/vV,  data  pres-jiiteu  by  Binding  ct 
al.  [15]. 


TABLE  3 


Material  parameters  for  Ml  at  20 °C  and  27.2°  C  using  (4.1)  and  (4.2) 


20  °C 

27.2°  C 

-;j,/2  =  a,  (gem  ') 

-3.0 

-1.0 

P  (gem  *) 

1.68 

0.54 

a2  (g cm'1) 

5.34 

1.77 

«I  +  «2 

2.34 

0.77 

n2/n, 

-0.11 

-0.12 
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5.  Second-order  correction  of  the  Trouton  viscosity 

We  turn  next  to  an  evaluation  of  the  coefficients  in  the  formula  (2.10) 
which  give  the  extcnsional  viscosity 

3[/i -r  (a,  *a3)i]  (5.1) 

at  low  rates  of  stretching.  Values  for  the  coefficient  or,  -r  a,  of  the  second- 
order  correction  of  the  Trouton  viscosity  are  given  in  Tables  2  and  3.  To 
complete  the  description  of  (5.1)  we  give  values  for  the  viscosity  \i  in  Table 
4.  The  elongational  stress  (2.10)  for  D1  at  25  °C  when  n2/nl  =  —0.10  is 

Tu  -7*22  ~  3(95i  -f  S3.0i2).  (5.2) 

Evaluation  of  the  same  stress  for  Ml  at  20  °C  gives 

Tj  i  —  Ta  ~  3(  29.5s  -f  2.34s2 ) .  (5.3) 

At  T—  21.2°  C  we  have 

Tj,  -  T22  ~  3(lls  +  0.11s2).  (5.4) 

Comparing  now  (5.2)  and  (5.3)  we  can  conclude  that  the  increase  of 


Fig.  10.  Comparison  of  cqn.  (5.3)  with  measurements  at  20  °C  of  Mikkclscn  et  al.  [7]  using 
opposing  jets. 


TABLE  4 


Values  of  the  viscosity  /t 


Fluid 

D1  at  25  °C 

9.5  Pa  s 

Ml  at  20°C 

2.95  Pa  s 

Ml  at  27.2°  C 

1.1  Pa  s 
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extensional  viscosity  at  second  order  is  more  than  50  times  greater  in  D1 
than  in  Ml. 

In  Fig.  10  we  have  compared  measurements  of  Mikkelsen  et  al.  [17],  using 
Fuller’s  method  of  opposing  jets,  with  our  expression  (5.3).  It  is  known  that 
the  opposing  jet  gives  a  stretch  value  for  the  extensional  viscosity  which  is 
about  1.4  times  the  Trouton  viscosity.  After  taking  into  account  the  expected 
shift  of  1.4  (see  Schunk  et  al.  [18]),  we  may  conclude  that  the  range  of  stretch 
rates  for  which  second-order  effects  are  negligible  is  the  same  for  the  exact 
formula  (5.3).  The  nonlinear  effects  in  the  experiments  undoubtedly  extend 
beyond  second  order  so  that  perhaps  there  is  satisfactory  agreement. 

6.  Wave  speeds,  effective  rigidity'  and  effective  relaxation  times 

The  shear-wave  speed  for  Ml  was  measured  with  the  wave-speed  meter  in 
our  laboratory.  The  meter  allows  us  to  determine  transit  times  for  an 
impulsively  generated  shear  wave  in  a  fluid  at  rest.  The  wave  traverses  the 
gap  between  two  concentric  cylinders.  Transit  time  measurements  must  be 
made  for  different  gap  sizes.  We  can  vary  the  gap  dimensions  in  our 
apparatus  between  0.25  and  12.8  mm.  A  shear-wave  speed  is  defined  if  we 
measure  one  and  the  same  transit  speed  over  a  range  of  small  gaps. 

The  theory  for  the  wave-speed  meter  is  given  in  Joseph  et  al.  [19].  The 
apparatus  and  the  measuring  technique  are  described  in  Joseph  et  al.  [20] 
and  in  detail  in  Riccius  [21], 

For  the  Ml  fluid  we  measured  transit  times  with  three  different  gap  sizes 
of  6.12,  9.38  and  12.2  mm.  The  transit  speeds  c  presented  in  Table  5  are 
determined  from  averages  over  eight  time  measurements  for  each  gap.  The 
shear-wave  speed  c  is  then  defined  as  the  average  over  the  three  transit 
speeds.  Examples  of  the  raw  data  for  each  gap  size  are  shown  in  Fig.  11.  The 
raw  data  arc  voltage  vs.  time  diagrams  which  are  recorded  on  an  oscillo¬ 
scope.  There  are  two  signals  for  each  gap  size;  each  one  is  adjusted  to  1.75  V 
before  the  measurement  is  started.  The  rapid  drop  of  the  first  signal 
corresponds  to  the  motion  of  the  outer  cylinder  which  is  caused  by  the 
impacting  kicking  device.  The  second  signal  represents  the  motion  of  the 

TABLE  5 

Wave  speed  c,  effective  rigidity  Gc  and  effective  relaxation  time  Ac  for  Ml 

Gap  c  T  c  Gc  \c  fi  p 

(10~3m)  (eras'1)  (°C)  (eras'1)  (Pa)  (10“3  s)  (Pas)  (kgm'3) 


6.12 

378.1  ±12.1 

23 

9.38 

403.4  ±33.6 

23 

12.2 

430.5  ±10.5 

23 

404.7 


14067  0.21 


3.0 


859 
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transit  times 


Fig.  11.  Evaluation  of  transit  times  from  raw  data  in  three  different  gap  sizes.  The  drop  in  the 
voltage  vs.  time  diagrams  shows  the  motion  of  outer  and  inner  cylinder  for  the  12.2  mm  gap 
(top)  (c  =  431  +  ll  cm  s"1),  for  the  9.38  mm  gap  (middle)  (c  =  403±34  cm  s"')  and  for  the 
6.12  mm  gap  (bottom)  (c  =  378±12  cm  s*1).  The  lengths  of  the  arrows  represent  the  lapse 
time  between  the  onset  of  motion  of  both  cylinders  in  each  gap  size. 


inner  cylinder  which  results  from  the  arriving  shear  wave.  How  the  transit 
time  is  evaluated  from  graphs  of  this  kind  is  indicated  through  the  arrows 
above  each  pair  of  signals.  The  lengths  of  the  arrows  represent  the  elapsed 
time  between  the  instances  of  onset  of  motion  of  the  outer  and  innei 
cylinder.  Figure  11  shows  the  development  of  the  transit  time  over  three  gap 
sixes.  From  here  and  from  Table  5  the  reader  can  verify  how  well  the  data 
satisfy  the  criterion  for  the  definition  of  a  shear-wave  speed  given  above. 

The  value  c  was  subsequently  used  to  evaluate  the  delayed  die-swell 
experiments  which  are  discussed  in  Section  7.  We  call  the  flow  “sub-”  or 
“supercritical”  whenever  the  local  velocity  is  smaller  or  larger  than  the  wave 
speed  c  that  was  measured  on  the  wave-speed  meter. 

In  order  to  make  a  comparison  with  other  fluids  tested  in  our  meter  we 
use  quantities  that  can  be  calculated  from  the  shear-wave  speed  c.  These  are 
the  effective  rigidity,  Gc  =  pc2,  where  p  is  the  density,  and  the  effective 
relaxation  time,  Xc  =  p/Gc,  where  p  is  the  zero-shear  viscosity  that  was 
determined  from  measurements  on  a  Rheometrics  System  Four  rheometer. 
We  believe  that  c,  Gc  and  Ac  are  relevant  for  the  dynamic  properties  of  flow 
and  for  steady  flows  which  change  type.  Gc  represents  slowly  decaying 
elastic  modes  that  characterize  wave  propagation  on  the  time  scale  of  the 
instrument  (0(lO-3-lO"4  s)).  Xc  is  an  important  estimate  for  the  decay 
time  of  these  modes. 

A  comparison  with  other  liquids  that  have  been  tested  on  the  wave-speed 
meter  is  especially  useful  (see  Riccius  [21]  for  a  complete  list  of  wave-speed 
data  obtained  until  February  1989).  The  signals  were  not  noisy,  facilitating 
unambiguous  evaluation  of  transit  times.  We  can  say  that  Ml  supports  wave 
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propagation  over  a  large  distance  at  high  wave  speeds  relative  to  most  other 
fluids.  In  comparison  with  fluids  of  similar  zero-shear  viscosity  the  large 
wave  speed  gives  rise  to  a  large  effective  rigidity  Gc  and  to  a  rather  short 
effective  relaxation  time  Xc.  The  relaxation  time  is  the  most  significant 
quantity  for  this  comparison.  It  effectively  reveals  whether  the  liquid  has  a 
relatively  elastic  or  Newtonian  response.  \c  plays  a  role  in  the  interpretation 
of  the  results  on  delayed  die  swell  (see  Section  7). 

7.  Delayed  die  swell 

Joseph  ci  al.  [22]  showed  that  there  is  a  critical  speed  of  extrusion  beyond 
which  the  swell  is  delayed;  instead  of  swelling  at  the  exit  lip,  as  it  does  in  the 
pre-critical  case,  the  extrudate  does  not  swell  at  the  exit  but  the  swell  begins 
somewhat  farther  downstream.  Many  photographs  are  shown  in  their  paper. 

A  sketch  of  delayed  die  swell  in  fluids  with  a  relatively  long  time  of 
relaxation,  say  Ac  >  10-2  s,  is  shown  in  Fig.  12. 

Joseph  et  al.  [22]  found  that  the  shape  of  the  jet  appears  to  correlate  with 
the  time  of  relaxation  based  on  our  wave  speed  which  was  discussed  in 
Section  6.  The  fluids  with  long  times  of  relaxation  are  more  elastic,  and  the 
delay  is  very  pronounced,  like  a  hydraulic  jump,  as  in  Fig.  12.  Fluids  with 
short  times  of  relaxation  are  more  Newtonian  with  larger  effective  viscosi¬ 
ties.  smoothing  shocks.  We  think  of  Jeffreys’  model  with  large  Newtonian 
viscosities,  large  retardation  times.  For  these  fluids  it  is  much  harder  to 
detect  the  delay.  The  delay  appears  as  a  change  of  curvature,  a  point  of 
inflection  in  the  jet  shape  which  is  introduced  at  criticality. 

In  Fig.  13  we  exhibit  a  sequence  of  photographs  which  shows  how  the 
shape  of  the  jet  of  Ml  changes  as  the  speed  of  extrusion  is  increased  past 
criticality.  The  extrusion  velocity  u  is  the  volume  flow  rate  divided  by  the 
area  of  the  pipe.  Ahrens  et  al.  [23]  have  shown  that  vorticity  of  the  Poiseuille 


Fig.  12.  Sketch  of  delayed  die  swell  in  a  solution  with  a  not-too-short  time  of  relaxation:  (a) 
subcritical,  (b)  supercritical.  The  speed  is  greater  than  the  wave  speed  before  the  swell,  less 
than  the  wave  speed  after  the  swell. 
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flow  of  an  upper-convected  Maxwell  model  changes  type  when  the  maxi¬ 
mum  speed  um  at  the  center  of  the  pipe  exceeds  the  wave  speed  c.  Since 
um  =  2  u  for  the  Maxwell  fluid,  theory  shows  that  a  hyperbolic  region 
appears  in  the  center  of  the  pipe  when 

M  =  um/c  =  2 u/c  >  1.  (7.1) 

The  procedures  used  to  obtain  the  values  of  u  and  c  are  exactly  the  same  as 
those  described  by  Joseph  et  al.  [22].  Photograph  13(a)  shows  a  subcritical 


Fig.  13.  Delayed  die  swell  of  Ml.  The  inside  diameter  of  the  pipe  is  3.175  mm.  (a)  u  =  312.88 
cm  s~\  M  =  2h/c  =  1.55,  subcritical.  (b)  u  =  357.58  cm  s-1,  M  =  1.77,  critical,  (c)  u  =  417.17 
cm  s-1,  M  =  2.06,  supercritical,  (d)  u  =  500.60  cm  s-1,  M  =  2.47,  supercritical. 


swell;  the  curvature  of  the  jet  is  one  signed.  Photograph  13(b)  shows  the 
condition  of  the  jet  slightly  after  criticality.  A  point  of  inflection  in  the 
curvature  appears  at  the  exit  lip.  We  can  see  this  on  the  videotapes  of  our 
experiments,  which  are  available  on  request,  as  a  sort  of  sudden  stiffening  of 
the  jet.  Photographs  13(c)  and  13(d)  show  the  jet  under  supercritical  condi¬ 
tions.  We  did  measurements  and  found  delays  in  pipes  with  different 
diameters. 

The  critical  conditions  for  delay  in  the  different  pipes  are  the  set  of  points 
belonging  to  Ml  shown  in  Fig.  14. 
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Fig.  14.  Mach  number  vs.  pipe  diameter:  o,  1.3%  CMC;  □,  9.8%  Elvacite;  o,  6%  PIB/D;  ©, 
2.5%  Polyox;  +,  12.1%  K-125;  a,  Ml. 


We  found  that  the  shapes  of  jets  of  Ml  were  essentially  identical  with 
shapes  in  other  fluids  with  short  times  X  =  j a/Gc  of  relaxation.  The  effective 
time  of  relaxation  for  Ml  is  0.00021  s.  We  refer  the  reader  to  Figs.  8.1  and 
8.2  of  the  paper  by  Joseph  et  al.  [22].  The  effective  time  of  relaxation  for 
PIBM,  1%,  and  PMMA,  1%,  are  listed  in  Table  1  of  Ref.  22  as  0.0020  s  and 
0.0059  s  respectively,  one  order  greater  than  the  2.6  X  10~4  s  for  Ml.  The 
jets  shown  in  Joseph  et  al.  [22]  were  extruded  downward,  in  the  direction  of 
gravity,  while  those  shown  in  Fig.  13  are  extruded  perpendicular  to  gravity. 
The  effect  of  gravity  in  the  vertical  jet  is  to  accelerate  the  jet,  to  stretch  it 
and  to  conceal  features  which  might  be  better  seen  in  a  horizontal  jet.  A 
PIBM,  1%,  jet  extruded  horizontally  would  probably  look  much  like  Ml. 

Ml  seems  to  be  close  to  a  Newtonian  fluid  with  small  non-Newtonian 
effects  which  are  often  associated  with  highly  viscous  liquids. 
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